Piezo-Magneto-Electric Effects in p-Doped Semiconductors 
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We predict the appearance of a uniform magnetization in strained three dimensional p-doped 
semiconductors with inversion symmetry breaking subject to an external electric field. We compute 
the magnetization response to the electric field as a function of the direction and magnitude of the 
applied strain. This effect could be used to manipulate the collective magnetic moment of hole 
mediated ferromagnetism of magnetically doped semiconductors. 



Antiferromagnetic dielectrics with inversion asymme- 
try exhibit the magnetoelectric (ME) effect, a phe- 
nomenon in which a static electric field induces a uni- 
form magnetization P,J3. Moreover, as first pointed 
out by Levitov et.al. l3!|, a kinematic magnetoelectric 
(kME) effect can also occur in ostensibly nonmagnetic 
conductors, with spin-orbit coupling, which lack a center 
of inversion symmetry. Unlike in dielectrics, in the case 
of conductors, the electric field induced magnetization 
density. Mi — aijEj, is necessarily accompanied by dis- 
sipation. Since M is odd under time reversal (T) and E 
is even, aij must be proportional to the relaxation time, 
a quantity related to the entropy production, making the 
process dissipative. In addition, since M is even under 
parity (P) while E is odd, is zero for parity invariant 
systems. The kME effect also vanishes in the absence 
of spin-orbit interaction. In 2D n-doped inversion layers 
with Rashba spin-orbit interaction, this effect has been 
predicted more than a decade ago 0, , but has been 
observed experimentally only very recently |^ Q ■ 

In the model used by Levitov et. al.^, the kME ef- 
fect originates from an electron scattering by impurities 
whose potential lacks inversion symmetry. As such, the 
effect is extrinsic and actually vanishes in the clean limit. 

In contrast, we present an analysis of hole-doped semi- 
conductor without inversion symmetry where the spin- 
orbit splitting of the p-band is intrinsic. In the absence 
of strain the system is both T and P invariant and hence 
no kME effect occurs. As we argue below, the shear strain 
induces a P-breaking term in the Hamiltonian which is 
responsible for the effect. 

There are several advantages to having a piezo- 
magneto-electric effect in 3D p-doped semiconductors 
(such as GaAs, GaSb, InSb, InGaAs, AlGaAS). Tech- 
nologically, engineering of different original strain archi- 
tectures is a common procedure in today's semiconduc- 
tor applications. By taking place in a 3D bulk sample, 
rather than in a 2D sample, this effect allows (with spe- 
cific strain configurations) full spatial manipulation of the 
magnetic moment. Most importantly, the effect occurs in 
p-doped semiconductors, and thus it allows manipulation 
of the direction of the collective ferromagnetic moment 
which develops in (p-doped) dilute magnetic semiconduc- 
tors. 

Within the spherical approximation 0,0], the effective 



Hamiltonian of a hole-doped semiconductor with spin- 
orbit coupling is described by the Luttinger-Kohn model 
in the spin-3/2 band: 
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where Si is the spin-3/2 (4x4 matrix) operator, 71 and 72 
are material-dependent Luttinger constants. The band 
structure consists of a doubly degenerate heavy hole band 
corresponding to k • S = ±3/2 and a doubly degenerate 
light hole band with k • S = ±1/2 (see inset of Fig.((Tl)). 
The above Hamiltonian is both P and T invariant. The 
strain, being a second order symmetric tensor , natu- 
rally couples to SiSj and to zero-th order modifies the 
original Hamiltonian Hlk by the term 

H^ = Dd{exx + f-yy + e-zz) + Due-ijSiSj, i,j^x,y,z (2) 

where and Z3„ are the usual hydrostatic and shear 
deformation potentials The modified Hamiltonian 

Hlk + H^ remains invariant under both P and T. Each of 
the two valence bands is still doubly degenerate. As seen 
in Fig. the strained Hamiltonian exhibits a finite 

energy gap between the heavy and light hole bands at 
zero momentum k=0. External electric field will cause 
a spin current 8], but no uniform magnetization. For 
semiconductors with inversion symmetry these are the 
only terms allowed at quadratic order in k. 

However, in the absence of an inversion symmetry cen- 
ter, the shear strain induces a P-breaking term linear in 
momentum 101: 



H' — XiSi, Xx 



C4 {^xyky 



(3) 



where Ay, Xz are obtained from Xx by cyclic permutation 
of indices and C4 is a material constant related to the 
interband-deformation potential for acoustic phonons. 
This term is responsible for the piezo-kME effect. Its 
origin can be traced back to the Kane's 8x8 model 
(2x2 for each the conduction and the split-off band, and 
4x4 for the valence band) (Fig. (^)), within which the va- 
lence band couples to both the conduction band and the 
split-off band. Upon straining, the zeroth order effect is 
the P-invariant coupling mentioned in the previous para- 
graph. At the first order, the conduction band \s) and 
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FIG. 1: Dispersion curves for InSb at 4 kbar stress on the 
[110] direction (tgy « 2 x 10~^) and k || [110] as measured 
by Seiler et al fll| . The strain splits the conduction and the 
valence bands. The Dresselhaus type inversion asymmetry 
is negligibly small on this scale. No splitting is observed for 
strain in the [001] direction. 



the valence bands \x), \y), \z) couple, and the matrix ele- 
ments between the valence and conduction band have the 
form exy{s\dxdy\z) (plus cyclic permutations) where |s) 
is the s-orbital and \z) is one of the p orbitals. Any other 
combination will not satisfy the selection rule. In sys- 
tems with inversion symmetry where the selection rules 
for L are satisfied, it is impossible to couple the spin-0 
(|s)) conduction states with spin-1 (ja;), \y), \z)) valence 
states through a spin-2 term (rank 2 tensor) (e^) and 
hence {s\dxdy\z) — 0. However, when inversion symme- 
try is broken, {s\dxdy\z) ^ as the L selection rule does 
not apply. We then obtain an 8 x 8 Kane matrix with 
the strain terms describing the interaction between va- 
lence and conduction bands. To find an effective 4x4 
hamiltonian for the valence band, one must project onto 
the valence band while taking into account the interac- 
tions with the conduction and the split-off band. The 
first term which appears in perturbation theory is the 
Hamiltonian JSJ. Reciprocally, a similar term will ap- 
pear in the conduction band, with the spin there being 
a spin-1/2 matrix. These terms have been observed ex- 
perimentally ^3 although recent evidence suggests other 
effects could also play a role 0. 

The piezo-kME effect can be easily understood as fol- 
lows: assume a material strained only along the [110] 
direction, such that exy ^ is the only nonvanishing 
shear strain component. Hence, the P-breaking term 



in the hamiltonian is H' = Ci{exykySx — £yxkxSy). 
This effectively corresponds to Zeeman coupling of hole- 
spins with a fictitious internal magnetic field Bx = 
Citxyky/ hb, By = -Ciexykx/^J-B, tJ-B being the Bohr 
magneton. Upon the application of an electric field 
along, say, the y axis, the average momenta become 
< kx >~ 0, ky K, eET/m where r is the momentum 
relaxation time. In turn, this gives < By >w 0, < 
Bx >~ CiexyeEr/mnB- The non-zero < Bx > field now 
couples to the spins and orients them along the x axis. 
This gives rise to a magnetization perpendicular to the 
electric field. Alternatively, the electric field along the x 
axis will induce magnetization along the y axis of equal 
modulus but of opposite sign to the previous one. This 
has recently been observed in the conduction band by 
Kato et.al. Ref.[^. Moreover, if we assume linear depen- 
dence on the relaxation time, r, and neglect the effect of 
parity conserving strain term, then the the form of is 
constrained by dimensional analysis alone 
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where n is the carrier density and the scaling func- 
tion ^ij{x,y) vanishes linearly with its first argument 
X. Based on the above argument, up to a sign, its com- 
ponents should be proportional to . 

We shall now justify the above claims. The static spin 
response to the d.c. electric field can be shown to be 
given by 
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where fis is the Bohr magneton and the retarded corre- 
lation function Q'^^^{uj) = Qp,i,{iuj lu + ir]), (77 ^ 0+), 
and 



dTe'^^{TS^{T)3,{0)). 



(6) 



For d.c. response only spatial averages of the spin and 
the current operators need to be considered above. The 
corresponding diagram is shown in Fig.Q. 

Since the strain splitting is typically small compared to 
the spin-orbit splitting at the Fermi surface (Fig^l, we 
can include its effects within (degenerate) perturbation 
theory. Utilizing the powerful mapping between spin- 
3/2 SU(2) and SO (5) representations pioneered by Mu- 
rakami, Nagaosa and Zhang ^8^, the unperturbed thermal 
Green's functions can be conveniently written as 

^ /I • \ 1 l-fsrf,(k)r,- 

2 -lujn + (1 + sA)e(k) 

where A = 272/71, e(k) — /2m, j = 1,...,5 and 
d(k) is a (spherical) unit vector in the 5-dimensional 
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FIG. 2: Top: Feynman diagram representing the kinematic 
Magneto-Electric efTect in 3D. Bottom: Kinetic equation for 
the vertex matrix in the ladder approximation. Here Vq'' is 
the velocity operator in the absence of impurities. 

space, equivalently d„i(k) = ^™2(^k,0k) where Y's 
are spherical harmonics and the angles are in k— space; 
Tj arc 5 Dirac gamma matrices 8j. Upon inclu- 



sion of spinless impurities with potential u(k), the full 
(impurity) Green's function G'(k, iw„) = G'o(k, -f 
I](k, zcjn)). Within the Born approximation the self- 
energy is E(k,iw„) = n.jnp /rfq \u{k - q)pGo(q, ici^„); 
nimp is the concentration of impurities. Finally, to lead- 
ing order in P breaking strain |l2l |. 

e(k, iujn) - [1 - G(k, iujn)H'ik)] G(k, iLUn). (8) 

Subsequently, all of the calculations will be carried out 
using ||SJ). 

As shown in Fig.|(2Jl, the finite frequency response func- 
tion © is given by 

Q^,{in) f dkTr [5,,n,(k, zc^, in)] , (9) 

where the trace is over the heavy /light hole spaces and 
where, as shown in Fig.||2Jl, the (4x4 matrix) vertex func- 
tion satisfies the kinetic equation (within the ladder 
approximation) 



U^Ckjiuj, ifl) = (/{k, iuj)Vuik)Q{k,iuj - ifl) + nmpQik, iuj) Jdq\u{k - q)fll^{q,iuj, in)Q{k, iu 
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(10) 



The velocity operator V^i(k) — dH{k) / dk^. Note that Q 
does not commute with li,^. As such we have 16 coupled 
integral equations to solve, one for each entry of the 4x4 
matrix. In the case of (S-function impurities u(k— q) = uq 
is a constant and the above integral equation is separable. 
Integrating both sides over k, it is easy to see that in the 
absence of parity breaking strain, the vertex correction 
vanishes [ij]. On the other hand, for finite strain, the 
vertex correction does not vanish, and we still have to 
solve a system of 16 coupled equations. 

However, to leading order in the strain it can be seen 
that all 16 equations decouple in the basis of the Clifford 
algebral Expanding the vertex matrix 

15 

n^(k, iuj, in) = (1^' " (ii) 

where the sum runs over all 16 elements [l4j and 
X"^ {k, iuj , iuj — in) is now an ordinary vector function. 
Since Tr[rA^iFB^i] is diagonal in A and B it is easy to 
see that 



where 



R^,{lLO,in) = UimpUQ I 2^ 
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and 



Vifiiuj, in) = \JdkTr [TAGik, iu;)V^{k)g{k, ito - in)] , 

(14) 

MA{ii^, in) = dkTr [TAQ[k, iuj)TAQ[k, iu - in)] . 

(15) 

With the known structure of the vertex matrix (|12|l . we 
can compute the response to the E^j^ field. Following the 
standard technique [15| we can perform the Matsubara 
summation, let in ^ n-\- irj, take the limit of — > and 
finally take the temperature T — s- to find 



^ rr r^v • riv ■ m c -ou '^^^ ^ ^ /dkTr[5^g-*(k)(K(k) + i?.)g"'^'^(k)], (16) 

-/ dkTr [rA5(k,«w)rBy(k, - zfi)] = M^(iw, 112)0^5 J 



Therefore, 

n^(k, iuj, in)^g{k, iu}){v^{k)+R^{iuj, in))g{k, iu-iniu) 



where the vertex matrix is given by the discontinuity 
of (|13|l . i?p = Rfi{iri, —irj). Finally, ignoring the inter- 
band transitions, (i.e. in multiple sums over s in Q we 
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keep only the same s), we get 
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where r is the momentum relaxation time, hb = 0.58 x 
10~^eV/G is the the Bohr magenton, n is the carrier 
concentration, and A = 272/71. Since 



(Ay, Az being obtained through cyclic permuations 
of x,y,z) are linear in the components k^^y^z, the factor 
dXi/dkj is a momentum independent, strain dependent 
tensor. For GaAs, C4 = 6*3/27? where C3 = 8 x lO^m/s 
|lO| is a measured constant related to the deformation 
potential for acoustic phonons while 77 = Ego /{Eg + 
Eso) = 0.183, Eg being the gap energy while Ego 
is the spin-orbit coupling energy for GaAs. For GaAs 
71 = 6.98,72 = 2.06, hence A = 0.59. This gives a 
value of Ca/H = 2.2 x lO^m/s. For n 



3 X lO^^cm- 



1%, in a generic sample of mobility /i = bOcrri^ /V -s 
the magnetization due (and perpendicular) to an elec- 



trie field is < S" >== 5 x W^E{V- 



). Under 



an electric field of £^ = W^V/rn the magnetization be- 
comes < S >= 6 X lO^^cm^'^. This corresponds to al- 
most < S >/n = 2% spin orientation efficiency. Since 
< S >^ kp ^ n^l'^ the spin orientation efficiency will 
grow for small hole concentration n. 

An interesting new application of the piezo-kME effect 
would be to manipulate the collective magnetization of 
the dilute magnetic semiconductors 01 ■ It is beheved, 
at least in the high mobility metallic regime, that the fer- 
romagnetism of Mn++ ions in GaMnAs is hole mediated. 
Within the k • p method [TtI , the coupling between 
the collective magnetization and the electric field induced 
spin polarization is of the order J « \eV [T^. H^. Thus, 
even if the total electric field induced magnetization rep- 
resents 1 Bohr magneton per 10'^ spins, the effective mag- 
netic field, felt by the Mn spins is tleff ~ ITTesla! 
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